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ABSTRACT 

' To unravel the formation mechanism and the evolutionary history of Elliptical Galax- 

PsJ I ies (EGs) is one of the goals of modern astrophysics. In a simplified picture of the issue, 

the question to be answered is whether they have formed by hierarchical merging of 
pre-existing sub-structures (maybe disc galaxies) made of stars and gas, each merg- 
^ \ ing event likely accompanied by strong star formation, or conversely, they originated 

from the early aggregation of lumps of gas turned into stars in the remote past via a 
burst-like episode ever since followed by quiescence so as to mimic a sort of monolithic 
process. Even if the two alternatives seem to oppose each other, actually they may 
■ concur to shaping the final properties of EG's as seen today. Are there distinct signa- 

' tures of the underlying dominant process in the observational data? To this aim we 

, have examined the line absorption indices on the Lick system of the normal, fiel d EGs 

of Tra gcr ( 1997) and the interacting EGs (pair- and shell-objects) of Longhetti et alJ 
l|200Clfl . The data show that both normal, field and interacting galaxies have the same 
scattered but smooth distribution in the H_a vs. [MgFe] plane even if the interacting 
O ' ones show a more pronounced tail toward high values. This may suggest that a 

\ common physical cause is at the origin of their distribution. There are two straightfor- 

C/3 . ward interpretations of increasing complexity: (1) EGs span true large ranges of ages 

and metallicities. The age youth is the signature of the aggregation mechanism, each 
event accompanied by metal enrichment. This simple scheme cannot, however, explain 
other spectro-photometric properties of EGs and has to be discarded. (2) The bulk 
population of stars is old but subsequent episodes of star formation scatter the EGs 
^ . in the diagnostic planes. However, this scheme would predict an outstanding clump at 

low values, contrary to what is observed. The model can be cured by supposing that 
the primary star formation activity lasted for a significant fraction of the Hubble time 
(5 Gyr^T^lS Gyr) accompanied by global metal enrichment. The "younger" galaxies 
are more metal-rich. The later burst of star formation should be small otherwise too 
many high H^ objects would be observed. Therefore, the distribution of normal, pair- 
and shell-galaxies in the H^ vs. [MgFe] plane is due to the global metal enrichment. 
Even though the above schemes provide a formal explanation, they seem to be too 
demanding because of the many ad hoc ingredients that have to be introduced. Fur- 
thermore they neglect the observationally grounded hint that the stellar content of 
EGs is likely enhanced in a-elements with [a/i^e] ranging from 0.1 to 0.4 dex. We pro- 
pose here a new scheme, in which the bulk dispersion of galaxies in the H^ vs. [MgFe] 
plane is caused by a different mean degree of enhancement. In this model, neither 
large age ranges nor universal enrichment law for the old component are required and 
the observed distribution along H^ is naturally recovered. Furthermore, later bursts 
of stellar activity are a rare event interesting only those galaxies with very high H^ 
(roughly > 2.5). Finally, simulations of the scatter in broad-band colors of EG's seem 
to confirm that the bulk stars have formed in the remote past, and that mergers and 
companion star formation in a recent past are not likely, unless the intensity of the 
secondary activity is very small. 

Key words: galaxies: elliptical - galaxies: evolution - galaxies: formation ~ galaxies: 
abundances 
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1 INTRODUCTION 

Determining the age of the bulk stellar content of Elliptical 
Galaxies (EGs) is basic to reconstructing the past history 
of star formation and to setting clues on the galaxy forma- 
tion process itself. Two competing scenarios are proposed: 
the monolithic picture in which a single dominant episode 
of star formation occurred in the past, ever since followed 
by quiescence and passive evolution (secondary episodes are 
however always possible); the hierarchical scheme, in which 
a series of mergers of smaller sub-units has taken place over 
the billions of years, each of those likely accompanied by star 
forming activity (mergers of inert objects with no additional 
star formation cannot be excluded). Current ideas concern- 
ing the two competing pi ctures of ga laxy formation hav e 
been critically reviewed bv lEUisI il998ft and iPeeble^ i2002l) . 
In short the present day view of structure formation suggests 
that roughly half of the large EGs were assembled at 2 < 1. 
However, observational evidence of variations with redshift 
of the color-magnitude, fundamental plane, size-magnitude, 
K-magnitude-redshift relations indicate that large EGs have 
already formed in the distant past at z ^ 2. 

In order to discriminate between the two schemes, the 
key question to be addressed and hopefully answered is 
what kind of signatures the stellar populations of a galaxy 
have inherited from the forming mechanism. The diagnos- 
tic to our disposal is mainly based on magnitudes, col- 
ors, integrated spectra, luminosity weighed line strength 
indices, estimates of chemical abundances and abundance 
ratios together with the spatial gradients in those quan- 
tities. As for the line absorption indices, gra dients in 
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have often been interpreted as indicating that the abun- 
dances of a-elements (Mg, O, etc.) with respect to iron are 
enhanced - [Q/-Fe]>0 - i n the central regions. Opposite 
concl usions however exist iKuntschnerl Il99a : iDavies et alJ 
I2OOII) . Furthermore, limited to the nuclear regions, indices 
vary passing from on e galaxy to another (iGonzaled Il993l : 
iTraeer et all l2000blal) . Looking at the correlation between 
Mg2 and (Fe) (or similar indices) for the galaxies in those 
samples, Mg2 increases faster that (Fe), which is once 
more interpreted as due to enhancement of a-elements in 
some g alaxies. In add i tion t o this, since the classical pa- 
per by iBurstein et alj lll988f) , the index Mg2 is known to 
increase with the velocity dispersion (and hence mass and 
luminosity) of the galaxy. Standing on this body of data 
the conviction arose that the degree of enhancement in 
a-element s ought to increase passing from dwarf to mas- 
sive EGs (iFaber et alJll992t lWorthev et aljll993 : iMatteuccil 
ll994ll997l:lMatteucci et alJll998h . 

Another point to consider is whether normal and inter- 
acting galaxies (i.e. shell- and pair-objects) show any sys- 
tematic difference in colors and line absorpt i on indi ce s. The 
subject was addressed by iLonghetti et al.l il998al lbl Il999l 
L2OOO), whose results will be shortly summarized below. 

Using this body of data, can we infer the age, the metal- 
licity and the degree of enhancement for the bulk stars in 
an EG? The task is not trivial because age and metallic- 
ity have similar effects on the spectrum of a galaxy, i.e. 



the spectrum and hence the colors of an old and metal- 
poor population may look like those of a young and metal- 
rich one, the so- called age m etallicity- degeneracy pointed 
out long ago bv iRenzini fc B uzzoni (1986). The degener- 
acy is also further complicated by effects both on the age 
and metallicity caused by the enhancement in a-elements. 
A promising way-out is perhaps offered by the system 
of line absorption indices introduced by the Lick group 

iWorthev ,1992; .Worthev et al. 199^ and reference therein) 

which have been extensively used to i nfer the age, metal- 
licity, and abundance ratios of EGs |B ressan et al j Il996l: 
Tantalo ct al. 1998; Kuntschnfi^^^^^^^998 ; Kunl^^^^ 
1998,; .J0r£cnsen 1999; Kuntschner 2000; Tragc r et alJ 
2000hllal: iKuntschner et aL 120011: i Vazdekis et al.. I2OOII; 
Davies et al.1 I2OOII: iPoggianti boOll: iMaraston et alJ 

20031: iTho mas et alj | 2003allbt iThomas fc MarastonI l200a 
Tantalo fc Chiosi 200j! 

Let us shortly summarize the key steps of those artic- 
ulated analysis, paying particular attention to the age for 
which often the conclusion that EGs show evidence of recent 
star formation either for the whole bulk stellar population 
or part it has been drawn. The age youth has often been 
taken as the signature of the hierarchical mechanism. 

(i ) Starting from the pioneering study of lGonzale'3 
il993h of "normal galaxies" (i.e. those with no sign of 
dynamical interaction, accretion, etc.) in the local Uni- 
verse, different groups analyzed the distribution of galax- 
ies in the v s. [MgFel diag nostic plane and oth ers of 
the same kind ("Buzzoni et al. 1992"; Buzzoni et al. 'l99j; 
Brcssan ct al. 1996; Tantalo 1998; Tantalo ct al. 1998) us- 
ing the Single Stellar Populations (SSP) approximation, i.e. 
the stellar content of a galaxy is reduced to SSP of suitable 
age, metallicity and degree of enhancement. The index 
is considered to be a good age indicator, whereas [MgFe] is 
considered to be most sensitive to the metallicity. In real- 
ity both indices are sensit ive to age, metallicity, and degree 
of en hancement (see also lGonzaleall993l : iTantalo fc Chiosil 
I2OO4I for similar remarks). Since EGs are more scattered in 
than in [MgFe] and do not follow the relation expected 
for c oeval old objects ma tching the Color-Magnitude rela- 
tion dBower et alJ Il992al lbl) . it was argued that II/3 traces 
the age of the bulk stars, thus suggesting that some galaxies 
are truly yo ung objects. Howeve r, a closer scrutiny of the 
problem led lBressan et all (|l_996') to suggest all galaxies in 
that sample should be old but with different histories of star 
formation. Some of them completed their stellar activity in 
the far past with no evidence of subsequent episodes. Others 
had a more prolonged star forming history, perhaps in re- 
current episodes of short duration. Looking at the galaxy to 
galaxy differences (5H/3 and (5 [MgFe] in the nuclear values and 
therr correlation with Eq, the suggestion was advanced that 
the global duration of the star forming activity gets longer 
at decreasing Eq (galaxy mass). This was later confirmed 
by dynami c al Tree-SPH simul ations of early type gala xies 
bv lKawatal lll999l l2001bllah andlChiosi fc Carrarol (1200211. In 
addition to this. iBressanet alJ (^96l ) and iLonriiBtti et alJ 
i200Cf) tried to understand the effect brought to the indices 
of a population of old stars by a recent episode of star forma- 
tion. The results is that indices like H/3 are strongly affected 
by even small percentages of young stars: as long as star 
formation is active they jump to very high values and when 
star formation is over they fall back to the original value 
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on a time-scale of about 1 Gyr. Other indices like Mg2 , (Fe) 
are much less affected even if the companion chemical en- 
richment may somewhat change them. In diagnostics planes 
like H/3 vs. (Fe) the galaxy performs and extended loop elon- 
gated towards the axis, thus causing an artificial disper- 
sion which could be interpreted as an age dispersion, whereas 
what we are really seeing is the transient phase associated to 
the temporary stellar activity. Consequently, the idea that 
the dispersion in measures the age of the last episode of 
star formation instead of the age of the bulk population was 
commonly accepted. Finally, EGs likely have mean metallic- 
ities in the range Zq ^ Z ^ 3Zq as ind icated by comparing 
theoretical m odels with data (see also iBressan et alj[l996l: 
laregriolll99^. 

(ii) iTantalo et alJ l)l998h examining the Gonzalez galax- 
ies, in order to infer from the simultaneous fit of If/j, Mg2, 
and (Fe) the age, the metallicity, and the enhancement factor 
[a/Fe], noticed that: (a) In a small group of galaxies, the nu- 
cleus appears to be more metal-rich, more a-enhanced, and 
younger (i.e. containing a significant fraction of young stars) 
than the external regions, (b) The galaxy to galaxy differ- 
ences SH13, 5Mg2, and 5(Fe) for the nuclear regions yield age 
differences for the last episode of star formation perhaps cor- 
relating with the galaxy mass. In other words, while there 
seems to be a sort of upper limit to the age of the bulk stars, 
traced by quiescent galaxies (no signs of rejuvenation), is 
seems also that the last episode of star formation took place 
closer and closer to the present at decreasing galaxy lumi- 
nosity (mass). 

(iii) Surprisingly, in the vs. [MgFe] plane, shell- 
, pair-galaxies (i.e. those with signs of dynamical inter- 
actions) and normal g alaxies share the same distribution 
iLonghetti et al.l I2OOOI) . There is, however, a group of pe- 
culiar galaxies with much stronger as compared to the 
normal ones. Does it mean that the scatter seen in this di- 
agram has a common origin, perhaps secondary episodes of 
star formation that can occur independently of whether or 
not a galaxy is interacting? This will be discussed in detail 
in Section 1^3] below. 

Thanks to these achievements, the vs. [MgFe] plane 
(and others of the same kind) for SSP of both solar and 
Q-enhanced abundance ratios soon became a popular tool 
to estimate the age, metallicity arid dey ee of enhancement 
of g alaxies (Kuntschnor & Davics 1998; Kuntschncr 199^; 
.Tc^rgen.sen 1999; Kuntschnor 2QQ0- Kuntschncr ct al. 200L: 
Vazdekis et al.l I2OOII: IPavies et al.l I2OOII: IPoeeianti et alJ 
2001f) . The advent of large libraries of SSPs with differ- 
ent chemical composi tions including the e ffect of differ- 
ent [a/Fe] ratios, e.g. ISalcisnich et alj (l2000l) . on one hand 
spurred the construction of corresponding libraries of line 
absorption indices, on the other hand made the many- 
indices fitting technique (three at least) of general use (e.g. 
[Trager et al. 2000b a; Maraston ct al. 2003; Thomas ct al, 
bnoaJhtlThomas fc Marastonll200.1lTantalo fc Chiosf20oll 
and introduced new techniques more sophisticated than the 
two-indices diagnostics. The solution for ages, metallicities 
and \a/Fe] ratios based on large samples of galaxies, e.g. 
the Trager (1997) list, once more yields large ranges of ages, 
metallicities, and abund ance ratios, as amply discussed by 
iTantalo fc Chiosil (|2Qq3). 

Although we have repeatedly mentioned that the large 
spread in age formally assigned to galaxies (EGs in partic- 



ular) on the base of the line absorption diagnostics has to 
be interpreted as referring to the last episode of star for- 
mation rather than the age of the bulk stellar population, 
quantifying the relative efficiency of the different star form- 
ing processes ha s still remained unsettled but for a few pre - 
liminary studies l)Bressan et alJl99fil : lLonghetti et al.l200(t l. 
The question is: how much a recent, minute episode of star 
formation, engaging a small fraction of the total mass, may 
alter the indices of an otherwise old population of stars in 
EGs? In this paper we reconsider the whole problem at the 
light of recent developments in the observational data and 
theory of population synthesis concerning the theoretical 
line absorption indices in which the effect of [a/Fe] ratios 
are taken into account. We intend in particular to investi- 
gate the effect of different histories of star formation and 
chemical enrichment from galaxy to galaxy and to explore 
in more detail the effects of late bursts of stellar activity. 

The paper is organized as follows. In Section |21 we 
shortly summarize the definition of line absorption in- 
dices and total enhancement parameter, mention the Fitting 
Functions in use, recall the main properties of stellar models 
and associated isochrones we have adopted to calculate SSPs 
and their integrated indices, and finally describe the key in- 
gredients of the library of stellar spectra we have used to 
calculate indices and broad-band colors of SSPs. In Section|^ 
we describe the indices for SSPs at varying metallicity, en- 
hancement parameter, and age. In Section |1| we present the 
technique to calculate the indices of composed objects made 
of at least two SSPs of different metallicity, enhancement, 
and age. In addition to this, we show the path in several 
two-indices-plane described by old galaxies rejuvenated by 
an episode of star formation of different intensity and age. 
In Section we present t he ob servational data we want to 
interpret, i.e. the iTrageil l)l997l) sample of normal galaxies 
complemented by the sample of iLonghetti et al . (200(|) for 
shell- and pair-galaxies. The analysis is based on the dis- 
tribution of galaxies in the diagnostic plane H/j vs. [MgFe]. 
Three different interpretations are explored: (i) The bulk 
stellar population of early-type galaxies span large ranges of 
metallicities and ages. This is equivalent to say that galax- 
ies can form over time interval comparable to the Hubble 
time, (ii) The bulk stellar population is very old (say 10- 
13 Gyr) in all galaxies but recent episodes of star forma- 
tions (bursts) engaging a small fraction of the total mass 
are present in most of them. Both explanations can hardly 
fit the distribution of galaxies along the H/3 direction unless 
ad hoc hypotheses are made for a sort of Universal Law of 
Metal Enrichment (ULME) . Since severe drawbacks of these 
explanations are either the large range of ages together with 
the ad hoc hypothesis of ULME and/or the otherwise un- 
avoidable existence of subsequent bursts of stellar activity 
in all galaxies, a third alternative is explored, (iii) The bulk 
population of stars in all galaxies is very old (once again 10- 
13 Gyr) but span large ranges of metallicity and global en- 
hancement factor. This is equivalent to say that each galaxy 
has its own history of star formation determining the chemi- 
cal properties. No ULME and no subsequent bursts of stellar 
activity are required to account for the observed distribution 
of the majority of galaxies. Only a few objects with unusu- 
ally strong H^ are likely in the bursting mode. In Sectional 
the last interpretation is critically discussed, the case of two 
template galaxies is examined in detail, and a few remarks 
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are made on the difFerent range of spanned by EGs in 
Coma and local vicinity. In Section |7| we address the ques- 
tion: how intense and old subsequent bursts of stellar activ- 
ity ought to be without violating the constraint imposed by 
the present-day mean broad-band colors of EGs? Finally, 
in Section |H] we summarize the results and presents some 
concluding remarks. 

2 THEORY OF LINE ABSORPTION INDICES 

2.1 Definition 

The technique to cal culate line-strength i n dices of SSPs 
is amp l y described in |Wor^]£3;^^yJl99j)j_Eressan_e t alJ 
Jl996l): iTantalo et al.l lll998l) : iMaraston et alJ J2003t) and 
iTantalo fc Chiosi J200j) . The reader is referred to those ar- 
ticles and references therein for a detailed description of the 
method. Here we limit ourselves to summarize a few basic 
points for the sake of clarity. Suffice here to recall that a line 
absorption index is constructed from the ratio Fi/Fc where 
Fi and Fc are the fluxes in the line and pseudo-continuum, 
respectively. The flux Fc is calculated by interpolating to the 
central wavelength of the absorption line, the fluxes in the 
midpoin ts of the red and blu e pseudo-continua bracketing 
the line dWorthev et alJll994l) . 

2.2 Fitting Functions 

Since the Lick s ystem of indices l|Burstein et al.1 11984 
iFaber et al]ll985l Tworthev et al.lll994) stands on a spec- 
tra library with flxed resolution of SA, whereas most of the 
synthetic spectra in use have a different resolution, to over- 
come the difficulty, the so-called Fitting Functions {J-J-) 
have been introduced. They express the indices measured 
on the observed spectra of a large number of stars with 
known gravity. Toff, and chem ical compo sition ([Fe/H]) as 
functions of these para meters ijW o^i gv et al. 1994) . In the 
following we adopt the lWorthevI lll99'^ TT. extended how- 
ev er to high temperature stars (Tcff ~10,000 K) as reported 
in lLonghetti et al.1 (Il9983). 

2.3 Integrated indices of SSPs 

They are derived in the following way. A SSP is described 
by an isochrone in the HRD whose elemental bins (fixed 
by AlogL/L© and A log Tcff) can be conceived as a "star" 
of suitable luminosity (gravity), Teff and metallicity. Each 
elemental bin is in turn populated by a number Ni of such 
stars given by 



the relative number of stars of the same type. We calculate 
the ratio 

where Ni is the above number of stars in the generic bin, 
and finally insert the ratio into the definition of each index 
as appropriate. 



2.4 a-enhanced chemical compositions 

Total enhancement. To distinguish the enhancement of 
individual species f rom the total one charac terizing a given 
chemical mixture, ITantalo fc Chiosil ll2004) introduce the 
parameter F defined as follows. Let us take a certain mixture 
of elements with mass abundances Xj of each species, and 
total metallicity Z (sum of the Xj for all elements heavier 
than He). The ratio of the mass abundance of an element 
with respect to Fe and to the Sun is given by 

AfeJ \AFe / \AFe J Q 

Let us now assume that the ratio [Xj /Xpe] of some elements 
is changed under the conditions that the abundance of Fe 
with respect to the Sun and the total metallicity Z remain 
constant. The mass-abundances of all species in the new 
mixture are accordingly scaled to a new value Xj . The total 
enhancement factor F is 

See ITantalo fc Chiosil ll2004h for ah other details. 

The iTripicco Belli (Il995h Response Functions. 

A method designed to include the effects of enhance- 
ment on line absorp tion indices has been suggested by 
iTripicco fc Belli lll995l TB95), who introduce the concept of 
Response Functions. In brief from model atmospheres and 
spectra for three stars of assigned effective temperature and 
gravity, i.e. a Cool-Dwarf (CD), a Turn-Off (TO), and a 
Cool-Giant (CG), they calculate the absolute indices Jq. 
Doubling the abundances Xi of the C, N, O, Mg, Fe, Ca, 
Na, Si, Cr, and Ti in steps of A[X,/H]=0.3dex they deter- 
mine the incremental ratios AIq/ A[Xi/ H] from which the 
Response Function Ro.3{i) for any index corresponding to a 
variation of the element t-th A[Xi/H] = -1-0.3 dex, can be 
derived 



(f){rn)dm 



(1) 



where nia and mt are the minimum and maximum star mass 
in the bin and 4>{m) is the initial mass function in number. 
We start from the typical star of a bin for which we may 
derive the indices from the TT. The indices are then in- 
verted to derive the fiux in the absorption line, F;*;, and 
in the pseudo-continuum F*^. Known the index for a single 
star, we weigh its contribution to the integrated value on 



Ro.3{i) = 



1 AJq 

Jo A[X^/H] 



0.3 



The Response Functions for the Cool-Dwarf, Turn-Off, and 
Cool-Giant stars constitute the milestones of the calibration. 

The mathematical algorithm for correcting an index 
from solar to Q-enh ance d element partitions is t aken from 
iTraeer et all (l2000b^ and ITantalo fc Chiosil J2004) to whom 
the reader should refer for all details. In brief, starting 
from the assumption that fractional variation of an index 
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Table 1. Chemical composition and [Fe/H] as function of F for 
the SSPs in use. 





F = 0. 


F = 0.3557 


F = 0.50 


Z Y X 


[Fe/H] 


[Fo/H] 


[Fe/H] 


0.008 0.248 0.7440 
0.019 0.273 0.7080 
0.040 0.320 0.6400 
0.070 0.338 0.5430 


-0.3972 
0.0000 
0.3672 
0.6824 


-0.7529 
-0.3557 
0.0115 
0.3267 


-0.8972 
-0.5000 
-0.1328 
0.1715 



to changes of the chemical parameters is the same as that 
for the reference index 7o, the following expression is derived 

^=^.{n[i+i^o.3W]*}-i (5) 

where Ro.^ij?) are the Response Function tabulated by TB95. 
Relation ^ can be applied under the obvious condition that 
I/Io > 0. 

The above algorithm is used to evaluate the frac- 
tional variations for the three calibrators, i.e. (A///)cd, 
{AI/I)to and (AI / I)cg- Since in general a star or elemen- 
tal bin along an isochrone will have effective temperature 
and gravity different from the ones of the three calibrators, 
and the fractional variations for these latter are also differ- 
ent, to evaluate the total fractional variation to be used for 
the particular star (isochrone bin) under examination, we 
linearly interpolate both in effective temperature and grav- 
ity among the fractional variations of the calibrators weigh- 
ing their contribution according to their distance from the 
current star in the HR-Diagram. For a det ailed discussion of 
this to pic and associated uncertainties see iTantalo fc Chiosi 
J2004l . 

2.5 Stellar models and isochrones 

The SSP indices are based on the Padova Library of stellar 
m odels and companion isochrones according to the version 
bv lGirardi et alJ J2000h and Girardi (2003, private commu- 
nication). This particular set of stellar models /isochr ones 
differs from the classical one by iBerteUi et ail Jl994) for 
the efficiency of convective overshooting and the prescrip- 
tion for the mass-loss rate along the Asymptotic Red Giant 
Branch (AGB) ph ase. For the reasons amply explained by 
iTantalo fc Chiosil 12004'). we prefer not to us e the more re- 
cent stellar models bv fSalasnich et all i2000ri in which the 
effect of a-enhancement is already included in the stellar 
opacity. Since indices are essentially a surface phenomenon 
in the sense that they are derived from the J^T linked to 
the stellar models only via the effective temperature, surface 
gravity, and iron content [Fc/H], details of the stellar models 
caused by patterns of abundances enhanced in a-elements 
are of minor relevance. The point has been made clear by the 
system atic analysis of the issue made bv ITantalo fc Chiosil 
^2004^■ Theref ore, the stellar models and companion SSPs 
bv lGirardi et a l. (2000) and Girardi (2003, private commu- 
nication) are fully adequate to the purposes of the present 
study. 

The stellar models extend from the ZAMS up to ei- 
ther the start of the thermally pulsing AGB phase (TP- 
AGB) or carbon ignition. No details on the stellar mod- 
els are given here; they can be found in iGirardi et alJ 



hood) and lGirardi et alj J2002l) . Suffice it to mention that: 
(i) in low mass stars passing from the tip of red giant 
branch (T-RGB) to the HB or clump, mass-loss by stel- 
lar winds is included according to the Rcimors (1975) rate 
with ?7=0.45; (ii) the whole TP- AGB phase is included in the 
isochrones with ages older than 0.1 G yr according to the al- 
go rithm of iGirardi fc Bertellil (^9^ and the mass-loss rate 
of IVassihadis fc Wood! Jl993^ ; (iii) four chemical composi- 
tions are considered as listed in Table Q 

2.6 Library of stellar spectra 

The l ibrary of stellar spectra is taken from IGirardi et alJ 
("2002^. ft covers a large range of the logTos - logg - and 
[M/ff] space. No details are given here. Suffice to mention 

• The basic spectra are from Kurucz A TLAS9 non- 
overs hooting models iCastelli et alJ Il997t iBessell et alJ 
Il998h complemented with: 

• Blackbody spectra for Toff>50,000 K; 

• iFluks et alj il994) empirical M-giant spectra, extended 
with synthetic ones in the IR and UV, and modified short- 
ward of 4000A so as to produce reasonable Tcs-{U — B) and 
Teff -(j3 — V) rela tions for cool giants; 

• lAUard et aP iioOO.'l DUSTY99 synthetic spectra for M, 
L and T dwarfs. 

The theoretical broad-band colors used in this study are 
in the Johnson-Cousins- Glass UBVRIJHK syste m, u sing fil- 
ter res ponse curves from lBessell fc Bretti lll98^ and lBessell 
Jl990l) . 



3 LINE ABSORPTION INDICES FOR SSPS 

In this study we essentially adopt the large grids of line 
absorption indices recently calculated bv .Tantalo fc Chiosi 
(2004) for SSPs over large ranges of chemical abundances, 
enhancement factors F, and ages. However we have slightly 
changed the ratio [Xei/Fc] for some specific elements. 
More precisely, the abundance ratios for solar-scaled and a- 
enhanced mixture s with F=0.35 and F=0.50 are taken from 
ITantalo fc Chiosil (|2004) but for Ti for which lower [Ti/Fe] 
are used. They are listed in Table (21 in which Columns (2), 
(3) and (6) show the abundance Aei of elements in logarith- 
mic scale, columns (4) and (7) list the ratio [Xei/Fe\, and 
finally columns (5) and (8) show the ratio [Xsi/HY . 

First of all and in view of the discussion below it is 
worth reminding the reader the definition of three indices 
that are commonly used but which do not belong to the 
original Lick system. They are (Fe), [MgFe] and [MgFe]' 

(Fe) = 0.5 X (Fe5270 + Fe5335) 



[MgFe] = \/Mgb X (0.5 x Fe5270 + 0.5 x Fe5335) 

^ The complete grids of stellar tracks, isochrones, magnitudes, 
are available on the web site http://pleiadi.pd.astro.it whereas 
SSPs colors, line absorption indices can be found on the web site 
http://dipastro.pd.astro.it/galadriel 
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Table 2. Abundance ratios for the solar-scaled and Q-enhanced mixtures adopted in this study. The abundances are the same as in 
[Tantalo & Chiosi 12004) but for Ti for which lower [Ti/Fe] ratios have been adopted: either [Ti/Fe]=0 (case a) or [Ti/Fe]=0.20 (case b). 
See the text for details). 





r = 


r = 0.35 


r = 0.50 




Element 


Aei 




I Fp J 






y fp J 


[^] 


Note 


O 


8.87 


9.37 


0.50 


HO. 1434 


9.57 


0.70 


4-0.2016 




Ne 


8.08 


8.37 


0.29 


-0.0666 


8.49 


0.41 


-0.0936 




Mg 


7.58 


7.98 


0.40 


(-0.0434 


8.14 


0.56 


+0.0610 




Si 


7.55 


7.85 


0.30 


-0.0566 


7.97 


0.42 


-0.0796 




S 


7.21 


7.54 


0.33 


-0.0266 


7.67 


0.46 


-0.0374 




Ca 


6.36 


6.86 


0.50 


HO. 1434 


7.06 


0.70 


-1-0.2016 




Ti 


5.02 


5.65 


0.00 


-0.3566 


5.89 


0.00 


-0.5013 


case a 


Ti 


5.02 


5.65 


0.20 


-0.1566 


5.89 


0.20 


-0.2201 


case b 


Ni 


6.25 


6.27 


0.02 


-0.3366 


6.28 


0.03 


-0.4731 




C 


8.55 


8.55 


0.00 


-0.3566 


8.55 


0.00 


-0.5013 




N 


7.97 


7.97 


0.00 


-0.3566 


7.97 


0.00 


-0.5013 




Na 


6.33 


6.33 


0.00 


-0.3566 


6.33 


0.00 


-0.5013 




Cr 


5.67 


5.67 


0.00 


-0.3566 


5.67 


0.00 


-0.5013 




Fe 


7.50 


7.50 


0.00 


-0.3566 


7.50 


0.00 


-0.5013 





and 



[MgFe]' = VMgb X (0.72 x Fe5270 + 0.28 x Fe5335) 

Going into a detailed description of the dependence of the 
line absorption indices for SSPs on age, chemical composi- 
tion, r, is beyond the scope of this study. The reader is re- 
ferred to[Tantalo & Chiosi (2004) who have discussed these 
matters in great detail. Suffice it to show in Fig.Qthe tempo- 
ral evolution of eight important indices, i.e. Mgb, Mg2, H/j, 
(Fe>, [MgFe], [MgFe]', NaD and C24668, for the following 
combinations of metallicity and F, namely Z=0.008 (solid 
lines) and Z=0.070 (broken lines), r=0. (heavy lines) and 
F = 0.35 (light lines). The age goes from 0.01 Gyr to 20 Gyr. 
It is worth noticing that the index Mg2 does not depend on 
F. 

Perhaps the most relevant result here is the variation of 
the indices with F. The easiest way of showing this is the 
vs. [MgFe] plane shown in Fig.|5| in which the range of values 
expected at varying metallicity (heavy solid lines), age (dot- 
ted lines), F and [Ti/Fe] (hatched areas) as displayed. Each 
hatched area is enclosed between the two SSPs with the low- 
est and highest metallicity in our sample, i.e. Z=0.008 (left) 
and Z=0.07 (right) and two lines of constant age (14 Gyr 
bottom and 2 Gyr top). Along each SSP four values of the 
age are marked, so that other lines of constant age can be 
drawn. Let us consider the hatched areas corresponding to 
r=0, 0.35 and 0.5, independently of the value for [Ti/Fe], 
either or 0.2. These cases will be discussed separately be- 
low. At increasing F both H/j and [MgFe] increase. While 
the effect on [MgFe] in somehow expected, the one on H/j 
is more difficult to understand. The explanation is in the 
strong R esponse Functions for elements like N, O, Mg etc 
found bv lTripicco fc BelJ ^1995^ and in the different way en- 
hancing those elements affects the ratio Fi/Fc- In presence of 
enhancement, the spectrum over the three pass-bands defin- 
ing the index H/jis more absorbed. However, absorption in 
the blue wing is larger than in the central band and in the 
red wing. Owing to this, the ratio Fi / Fc gets smaller so that 
H;3 gets stronger. This trend is also confirmed by the new 




5 10 15 20 5 10 15 20 

T(Gyr) T(Gyr) 



Figure 1. Evolution of eight indices (Mgi,, Mg2, H^, (Fe), 
[MgFe], [MgFe]', NaD and C24668) as function of the age. The 
heavy lines show the indices for solar-scaled partition of elements 
(r=0), whereas the thin line s show the same but f or the parti- 
tion of a-enhanced taken from lSalasnich et al.H200(f) (r=0.3557). 
Only two metallicities are displayed for the sake of clarity, i.e. 
Z=0.008 (solid lines) and Z=0.07 (dashed lines). 

Response Functions calculated bv lTantalo et all (|20o3) us- 
ing high resolution spectra (lA). 

There is another important poi nt to be addressed , which 
has already been touched upon bv lTantalo fc Chiosil (|20o3) 
and it is shown in Fig. |5| by the cases with the same F, 
either 0.35 or 0.5, and different [Ti/Fe]. In brief, at given 
metallicity Z, total enhancement factor F, and list of en- 
hanced/depressed elements, the same F c an be obtained by 
many patterns of [Xei/Fe]. For instance iTantalo fc Chiosil 
i2004l) have shown that at given F indices like are very 
sensitive to the abundance ratio [Ti/Fe]. In their calculations 
several values of [Ti/Fe] have been explored, i.e. [Ti/Fe]=0, 
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g [Ti/Fe]^0,Z & r-0,50 - 
^ [Ti/Fe]=0.0 & r=0.50_ 




0.070 



1 2 3 4 5 6 

[MgFe] 

Figure 2. The vs. [MgFe] plane for different combinations 
of r=0, metallicity, and [Ti/Fe] as indicated. Eacli hatciied area 
is enclosed between the two SSPs with the lowest and highest 
metallicity in our sample, i.e. Z=0.008 (left) and Z=0.07 (right). 
Along each SSP four values of the age are marked so that the lines 
of constant age can be drawn. Four values of the age are, i.e. 2, 3, 
8, and 14Gyr. The different hatched areas correspond to different 
combinations of F and [Ti/Fe] as indicated in the legend. Looking 
at this diagram one can easily single out the separate effect of the 
four parameters: age, metallicity, F and [Jf /Fe] ([Ti/Fe] in this 
case). See the text for more details. 

0.20 and 0.63. The first choice implies that Ti is not consid- 
ered as an Q-eleme nt, the second one is based on the mean 
value measured bv lGratton et all (|20o3) for galactic metal- 
poor stars with accurate parall axes, wherea s the l ast choice 
comes from the old estimate bv lRvan et al.l il99lll for stars 
of the same type. In Fig. |5] we show the effect of varying 
[Ti/Fe] from to 0.2. Test calculations show that other el- 
ements like O, Mg, Ne, Ca etc. have a similar effect. The 
explanation is the same as before. 

Since there is no unique pattern of abundance ratios as 
clearly shown by the above mentioned observational data, 
at given Z and F a sort of natural width for the indices has 
to be expected. The natural width we are talking about is 
provided by the difference [lyx^jFe] - hx^i/F^y) age,r,z- 
is soon evident that models with the same Z, V and age 
may have significantly difi'erent and [MgFe] (indices in 
general) depending of the detailed pattern of [Xei /Fe] in use. 



4 COMPOSING SSPS 

Since the stellar content of a galaxy can be approximated 
by a manifold of SSPs of different age and metallicity, each 
of which weighed on the star formation rate, the indices of 
composite stellar populations can be derived in simple fash- 
ion without making use of sophisticated population synthe- 



sis techniq ues requiring chemo- spectro-photometric models 
of galaxies jTantalo et alJll998l) . 

Given the generic index /f ^^(F, Z, F, t), with the aid of 
eqn. Q we derive the integrated fluxes of the SSP in the 
absorption feature and pseudo-continuum 

Tl{Y,Z,V,t)=Y.K^N^ (6) 

i 

and 

Tc{Y,Z,V,t) = Y.K^Ni (7) 

i 

Let us consider a composite object made of several SSPs 
with assigned composition, total enhancement and age - 
{Yj , Zj ,Tj ,tj) - each of which intervening in suitable per- 
centages /3j. This parameter measures the intensity of star 
formation assigned to each component, in other words it is 
the fraction of the galaxy mass engaged in each SSP. The 
integrated fluxes are 

j-comp ^ ^ j-^ ^y. ^ Zj ,Vj,tj) (8) 

j 

and 

^ = ^/3,^.,,(F„Z„F„t,) (9) 
i 

The total mass is the sum of the SSP mass over all the com- 
ponents and re-normalization of the fluxes can be applied if 
required. The composite indices are then derived by insert- 
ing the above integrated fluxes into eqn. Q. 

For the sake of illustration we consider the case of an 
old galaxy, represented by a SSP with solar composition and 
no enhancement of a-elements (Z=0.019 and F=0), which 
at the age of 10 Gyr suffers a burst of stellar activity rep- 
resented by another SSP of the same Z and F. The old 
galaxy is labelled by j—1 and the burst component by j=2. 
Three cases of the burst intensity are examined, i.e. /3i=0.98 
and /32=0.02, /3i=0.90 and /32=0.10, /3i=0.80 and /32=0.20. 
Once the burst has occurred the photometric properties of 
the composite object are described with the age step of the 
young SSP in order not to loose in the time resolution of 
magnitudes, colors, and indices. In simulations of this type, 
both the old and the young component are let evolve with 
time from the burst epoch to the present. The results are 
shown in the panels of Fig. |21 and Fig. |1| The left panel 
of Fig. 13 is the color (B-V) vs. age (in Gyr); whereas the 
right panel is the vs. age. The left panel of Fig.^Jis the 
plane (U~B) vs. (B-V), whereas the right panel is the vs. 
[MgFe] plane. For the sake of clarity, the burst is displayed 
only for ages older than 0.1 Gyr because its path in the two 
colors and/or two indices planes for younger ages is quite 
complicated and difficult to describe in a simple fashion. In 
the left and right panels of Fig. |21 the age structure of the 
bursting mode is clear. The color (B-V) and index are 
suddenly rejuvenated at the onset of the burst: (i) (B-V) 
may get very blue, the peak value depends on (32, and then 
fade down toward the typical color of an old object as the 
age increase from 10 Gyr up to 16 Gyr. The same for H^. 
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Age{Gyr) Age{Gyr) (B-V) [MgFe] 



Figure 3. A burst of star formation of different intensity (as 
indicated) is let occur in an ideal galaxy represented by a SSP 
with Z=0.019 and r=0 at the age of 10 Gyr. The composite object 
is then followed up to the age of 16 Gyr. Left Panel: (B-V) vs. 
age. Right Panel: vs. age. 



(ii) In all simulations both and (B-V) keep memory of 
the bursting activity for long time. Both indeed tend to be 
lower and higher, respectively, than the pure passive evolu- 
tion indicated by the solid heavy line. In the planes (U-B) 
vs. (B-V) and H;3 vs. [MgFe] of Fig. ^ the composite ob- 
ject performs wide loops, the extension and width of which 
are functions of the intensity 132. The complete path of the 
bursting object is better shown in Fig. |3 below. 

Another interesting plane to look at is the superposi- 
tion of two SSPs, one old and the other young (in different 
percentages as above, but in which the age of the old SSP 
is kept fixed whereas that of the young one is let vary from 
0.01 Gyr to the age of the old component in suitable steps. 
These simulations allow us to span the whole range of pos- 
sible combinations of ages for the two SSPs and the whole 
range of values in any two-indices plane at varying intensity 
and age of the burst. In other words, any point in these plane 
is a picture of the composite galaxy taken at a certain age 
of the burst. Simulations of this type are shown in Fig. |K| 
for the typical age of the old SSP of 13 Gyr. The appar- 
ently strange behavior of the SSP path at increasing age of 
the burst deserves some explanation. The bell-shaped tren d 
of has already been explained bv lBuzzoni et all il994) . 
In brief for strong absorption (as in young SSPs and/or 
burst of star formation), the Stark wings of the feature over- 
flow the Lick wavelength window, and enter the side bands 
depressing the pseudo-continuum. As a result, the index 
peaks when the SSP is dominated by B5-A0 stars and fades 
for earlier and later spectral types. This is mirrored in the 
path of the composite SSP with the contribution from each 
component weighed on its percentage. It is worth noticing 



Figure 4. The same as in Fig. |2l Left panel: the (U— B) vs. 

(B-V) plane. Right Panel: the vs. [MgFe] plane. 




Figure 5. Composite SSPs: young SSPs of any age and metallic- 
ity Z=0.040 are superposed to an old SSP with the age of 13 Gyr 
and the metallicity Z=0.019. The age of the old SSP is kept con- 
stant. Different percentages {132) for the young component are 
considered as indicated. The percentage of the old component 
is /9i=l — /32- The various panels show four different diagnostic 
planes. 
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that a galaxy caught at the very early stages of its burst- 
ing activity would appear as an object with unusually low 
indices. We will come back to this later on. 



5 INDEX-INDEX DIAGNOSTICS: AGE OR 
CHEMISTRY? 

The sample of galaxies we intend to analyze is derived from 
two sources o f data, namely the "/DS Pnstine" catalog by 
iTrae er' f 1997*) for the central regions of normal galaxies - it 
contains also th e Gonz alez (1993) list and the catalog by 
iLonghetti et al.l ll200Cf) for pair- and shell-galaxies (i.e. only 
objects with clear signs of interaction) . The latter sample is 
given in Table|3f . The H/j vs.[MgFe] plane for all the galax- 
ies in question is shown in the left panels of both Figs.|S]and 
Fig.|7|to facilitate the comparison with the theoretical simu- 
lations we are going to present. It is worth calling attention 
on the striking similarity between normal and interacting 
galaxies, and the very smooth but steeper slope of the data 
as compared to that of SSPs of given metallicity (see Fig. 
I^J. In the following we will refer to it as the nearly vertical 
distribution along the axis. 

5.1 A toy galaxy model 

iBressan et alJ (liggef) and ' Longhetti et al.l J2000 h suggested 
that the distribution of galaxies in the vs. [MgFe] plane 
may reflect secondary episodes of star formation superposed 
to an old stellar component. In the following, we intend 
to explore further thi s idea by means of Mo nteCarlo sim- 
ulations based on the ILonghetti et al.l i2000f) toy model of 
galaxy evolution. 

The complex star formation history of an early-type 
galaxy is reduced to a burst of relatively recent star forma- 
tion superposed to the bulk population made of old stars. 
These latter are in turn represented by a SSP whose age is 
randomly selected between the ages Ti and T2. The young 
stellar component, formed during the recent burst of star 
formation, is represented by a SSP whose age is randomly 
selected between T2 and a lower limit T3. Typical values 
for the age limits are: Ti=13, T2=10, and r3=0.1 Gyr. The 
time-scale of the two star forming event are always assumed 
to be short compared to all other relevant timescales (ages, 
and Hubble time). 

Since we are interested in guessing the minimum thresh- 
old above which the secondary episode gets importance in 
affecting the line strength indices, we will consider only the 
case in which the secondary episode involves a minor frac- 
tion of the galaxy mass. Typical values for the strength of 
the primary and secondary burst /3i and f32, respectively, 
measuring the percentage of the total galaxy mass turned 
into stars, are /3i=0.98 and /32=0.02. 

The metallicity Zi of the old stellar component is ran- 
domly selected between 0.008 (50% of the solar value) and 
0.07 (3.25 the solar value). We have also built a set of sim- 
ulations in which the metallicity of the old component is 

^ We remind the reader that the values of [MgFe] listed in Table 4 
of^onghetti et al. (2000) have been calculated using the wrong 
expression [MgFe]= y^Mgb X Fe5270 + Fe5335/2 



forced to linearly increase from 0.008 and 0.07 over the age 
range Ti to T2 . The metallicity Z2 of the young component 
is randomly chosen over the whole range (i.e. from Z=0.008 
to Z—0.07), thus simulating the widest range of possibili- 
ties, going from acquisition of external less processed gas to 
chemical enrichment during the burst. 

The simulations are first performed for solar partitions 
of elements, i.e. F^O and then for F 7^0 (another dimension 
is added to the problem). 

Finally, random errors are applied to the model in- 
dices to bette r simulate the observations. Using the data by 
iTrageil il997l) we calculate the mean relative errors {A///)o 
as reported in Table 3] The error affecting an index is ran- 
domly evaluated according to 

where e is a random number between and 1. 



5.2 Large scatter in age and metallicity of the 
bulk population 

This is the simplest interpretation of the distribution of the 
galaxies in the vs. [MgFe] plane. Neglecting important 
effects due to enhancement in a-elements and recent stel- 
lar activity, matching the observational range of the data 
would require the bulk population of early type galaxies be- 
ing formed in different epochs from galaxy to galaxy over 
a time-scale comparable to the Hubble time. In the cen- 
tral panel of Fig. |S| we show a simulated sample of 100 
objects for which the following parameters are assumed: 
ri = 13Gyr, r2=2Gyr, /3i = l, no increase of the metallic- 
ity in this time interval, no later bursts of stellar activity, 
and finally r=0. The models galaxies distribute along the 
SSP lines of different metallicity according to their age. The 
simulation significantly differs from the observational data 
shown in left panel of Fig.|H] Even if this view could be fitted 
into the classical hierarchical scheme of galaxy formation, it 
can be hardly sustained because it would predict spectro- 
photometric properties not fully compatible with the obser- 
vational data for EGs. 



5.3 Random bursts of star formation 

Bursts of star formation (from one to several) superposed 
to an old dominant population of stars seem to be more 
plausible and yield better results. Galaxies are conceived as 
old, nearly coeval systems, their population being approxi- 
mated by a single SSP with age between 10 and 13 Gyr; this 
age range agre es with the current age estimate of EGs in 
rich clusters |Bower et al, 19981). A burst of stellar activity 
is added at an age randomly chosen in the interval T2-T3. 
Two different prescriptions for the metallicity are adopted 
as described in Section lS.ll above. Finally all the simulations 
are for r=0. The maximum intensity of the superposed burst 
amounts to 2% of the total mass, i.e. /3i=0.98 and /32=0.02. 

The simulations of this type are made at increasing 
complexi ty. Since they essentia lly confirm what already 
found bv ILonghetti et all (I2OOOI) we limit ourselves to dis- 
cuss the results and to highlight the point of disagreement 
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Figure 6. Observational data and MonteCarlo simulations in the vs. [MgFe1 plane. L eft Panel: The normal galaxies bv lTrageJ ll997f) 
indicated by the open triangles and the shell- and pair-galaxies of lLonghetti et all i200(f) indicated by the filled squares. Note the smooth 
and broad distribution of both groups of galaxy. Central Panel: MonteCarlo simulations for the set of parameters Ti=13 Gyr, T2=2 Gyr, 
/3l=l, no additional bursts of star formation. The metallicity is randomly chosen in the interval 0.008^ZsC0.07, no enhancement of a- 
element is considered (r=0). The solid lines are two SSPs with Z=0.008 (left) and Z=0.07 (right). Four values of the age, i.e. 2, 3, 8 
and 13 Gyr are marked. The thin dotted lines are the loci of constant age. The very large age range is required to get full coverage of the 
range in spanned by the data. Right Panel: MonteCarlo simulations according to the recipe: ri=13Gyr, r2=5Gyr, random burst 
of star formation activity at the time 0.1^ T3 ^10, the same interval for the minimum and maximum metallicity and the same V as in 
the previous simulations. The percentage of the young population is /32 =0.02. The bulk population formed in the time interval T\ to T2 
obeys the universal law of metal enrichment. See the text for details. All other symbols have the same meaning as in the Central Panel. 



Table 3. Basic data for the iLonghetti et al.l l200Ci) sample of shell- and pair-galaxies. The velocity dispersion Eq is km/s. 









Pair-Gal 


axies 














Shell-Gal 


axies 








Name 




Mg2 


Mgb 


Fo52 


Fo53 


[MgFo] 


So 


Name 




Mg2 


Mgb 


Pe52 


Fe53 


[MgFe] 


So 


RR24a 


-2.60 


0.17 


2.49 


2.22 


2.32 


2.38 


101 


N813 


2.29 


0.25 


4.19 


3.33 


2.28 


3.43 


208 


RR24b 


-17.37 


0.09 


2.01 


0.74 


1.10 


1.36 


92 


N1210 


1.44 


0.26 


3.94 


2.92 


2.46 


3.26 


180 


RR62a 


1.15 


0.13 


2.30 


1.57 


1.82 


1.97 


81 


N1316 


2.01 


0.26 


4.02 


2.84 


2.62 


3.31 


250 


RRlOla 


1.86 


0.23 


3.53 


3.10 


2.76 


3.22 


162 


N1549 


1.71 


0.29 


4.46 


2.85 


2.64 


3.50 


225 


RRlOlb 


2.35 


0.23 


3.67 


2.83 


2.35 


3.08 


210 


N1553 


1.40 


0.30 


4.45 


3.05 


2.79 


3.60 


188 


RRlOSa 


1.32 


0.27 


4.23 


3.21 


2.26 


3.40 


130 


N1571 


1.66 


0.28 


4.31 


3.24 


2.48 


3.51 


224 


RR187b 


2.87 


0.22 


3.40 


2.69 


2.34 


2.92 


144 


N2865 


3.12 


0.22 


3.28 


2.34 


2.22 


2.22 


208 


RR210a 


0.80 


0.33 


4.97 


2.71 


2.47 


3.59 


281 


N2945 


0.38 


0.26 


4.80 


2.62 


2.21 


3.40 


212 


RR210b 


1.50 


0.31 


4.62 


2.99 


2.44 


3.54 


212 


N3051 


1.15 


0.30 


4.82 


2.74 


2.45 


3.54 


212 


RR225a 


1.76 


0.34 


4.96 


3.41 


2.40 


3.80 


357 


N5018 


2.68 


0.22 


3.31 


2.89 


2.24 


2.91 


247 


RR225b 


1.18 


0.31 


4.42 


2.80 


2.59 


3.45 


186 


N6776 


1.92 


0.25 


4.39 


2.67 


3.15 


3.57 


243 


RR278a 


-1.67 


0.20 


4.14 


2.20 


2.44 


3.10 


184 


N6849 


1.32 


0.28 


4.34 


3.03 


2.67 


3.52 


198 


RR282b 


1.53 


0.33 


5.17 


2.78 


2.16 


3.57 


290 


N6958 


1.66 


0.26 


4.08 


2.80 


2.29 


3.22 


223 


RR287a 


1.61 


0.28 


4.57 


3.08 


2.43 


3.55 


205 


E1070040 


2.30 


0.25 


3.97 


2.98 


2.74 


3.37 


155 


RR297a 


1.97 


0.27 


4.08 


2.89 


2.32 


3.26 


177 


E3420390 


1.32 


0.33 


4.77 


3.15 


2.75 


3.75 


327 


RR297b 


1.17 


0.26 


4.81 


2.92 


2.60 


3.64 


147 


N7135 


-0.41 


0.31 


5.51 


2.81 


2.38 


2.38 


191 
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91 


E5380100 


1.70 


0.31 


4.20 


3.01 


2.10 


3.28 


178 


RR381a 


1.64 


0.31 


5.01 


2.63 


2.62 


3.63 


276 


















RR387a 


2.20 


0.29 


4.48 


3.08 


2.74 


3.61 


187 


















RR387b 


1.83 


0.31 


4.51 


3.13 


2.55 


3.58 


215 


















RR397b 


1.33 


0.29 


4.39 


3.21 


2.87 


3.65 


198 


















RR405a 


1.45 


0.28 


4.68 


2.92 


2.39 


3.52 


192 


















RR405b 


1.89 


0.34 


4.81 


2.72 


2.76 


3.63 


207 


















RR409a 


1.89 


0.29 


4.51 


3.02 


3.28 


3.77 


187 


















RR409b 


1.17 


0.31 


4.56 


2.89 


2.36 


3.46 


150 



















Star Formation History in Early- Type Galaxies. I. The Line Absorption Indices Diagnostics 11 



with the observational data without showing any simulation 
in detail: 

(i) Stronger bursts of star formation (i.e. engaging more 
than 2% of the mass) are nor suited as they would predict 
too high values of and too many young objects. 

(ii) The expected distribution of objects with respect 
to the H/3 index is at variance with the observational one. 
Indeed models of this type predict a bimodal distribution, 
whereby the old galaxies (those for which the burst is al- 
most as old as the bulk of their stellar populations) clump 
together in the lower portion of the diagram, whereas the 
"young" objects (those with very young bursts) form a tail 
extending to high values of H/j. 

(iii) Nevertheless, the burst alone cannot explain the 
smooth distribution observed at low H/3 values. This is be- 
cause the index of a stellar population for which the 2% 
of the mass is composed by "young" stars and the remaining 
98% by an old component, reaches the observed high values, 
but, fading very rapidly with the age of the "young" compo- 
nent, has low probability to match the intermediate observed 
values. Slowing down the index decrease (corresponding to 
the aging of the "young" component) by increasing the per- 
centage of mass involved by the young burst produces a 
uncomfortably large fraction of objects in the upper part 
of the diagram. A complex interplay between burst inten- 
sity and mean age of the stellar population should then take 
place, with the old bursts being on average stronger than 
the recent ones. It must be said, however, that when the 
burst itself is larger than a few percent of the total mass, 
the definition of the average age of the bulk of the stellar 
population becomes a problem. 

(iv) The difficulty is partially cured assuming that the 
old population has an average age spreading over a signifi- 
cant fraction of the Hubble time (say down to T2 ~5Gyr). 
This is meant to indicate that either the object has been 
growing for such a long time with a low star formation rate, 
or that its major star formation activity was not confined 
to an early epoch. The young component is left to occur. 
It appears immediately that the observed smooth distribu- 
tion in the H^ index together with the young tail can be 
much better reproduc ed with such a kind of models, see 
iLonghetti et al.l hOOd) for details. 

(v) Another problem arises if the metallicity is ran- 
domly selected. The distribution of the models galaxies in 
the left panel of Fig. |S| strictly follows the path of a SSP, 
whereas the data run much steeper. We take this point to 
invoke the existence of a relation between the age of the bulk 
stellar population and its average metallicity and to apply 
the ULME. The simulations better fit the data. 

The right panel of Fig. |H| shows our final experiment 
incorporating all the hints we have been discussing so far. 
The simulations are based on the following parameters and 
assumptions: Ti^lSGyr, r2=5Gyr, r3=0.1Gyr, /3i=0.98, 
/32=0.02, average metallicity of the bulk of the stellar popu- 
lation forced to linearly increase from Z=0.008 to Z=0.070 
over the time interval Ti to T2. Thanks to the combined 
effect of the large age and metallicity spread for the bulk 
population, the distribution in H^ vs. [MgFe] plane is nearly 
vertical. The larger range of metallicity is necessary to main- 
tain a significant dispersion in [MgFe], because age differ- 
ences tend to compensate metallicity differences. If the latter 
interpretation is correct, it suggests that young early- type 



galaxies in the field are on average more metal-rich than 
old systems, with an average metallicity gradient of about 
Alog(Z)/Alog(t)~ —0.7, the latter value being very depen- 
dent on the younger limit of the bulk age. 

The main conclusion out of these simulations is that in 
addition to the mass dominating old population, which how- 
ever has to be built up over a large time interval and under 
a suitable age-metallicity relationship, sprinkles of stellar 
activity in the recent or very recent past ought to consid- 
ered in nearly all galaxies to reconcile theory and observa- 
tions. However, that all galaxies have to go through recent 
star forming activity is perhaps too demanding and other 
alternatives should be explored. 

5.4 The Q-enhancement alternative 

The resuhs obtained bv iTantalo fc Chiosi J2004 for SSPs 
of the same age and metallicity but different degrees of en- 
hancement in a-elements offer a third plausible explanation. 
The bottom line of the model is best explained by comparing 
the theoretical models in the H/3 vs. [MgFe] plane of Fig. |S| 
(central and right panels), with the observational data (left 
panel). Several points are soon evident: 

(i) The majority of galaxies (those with H/3^2) are fully 
compatible with being very old objects of the same age (say 
about 13 Gyr) but a different degree of enhancement in q- 
elements going from F^O to r=0.5 (taking the so-called nat- 
ural width caused by a possible variation in single elemental 
species into account). As a matter of fact, an old galaxy 
(say a 10-13 Gyr object) being shifted to higher H/3 by high 
F and/or [Xej/Fe] ([Ti/Fe] as a prototype) could lie in the 
same region occupied by a galaxy of significantly younger 
age and solar abundance ratios. At least part of the scat- 
ter along the H^ axis could be due to a different degree of 
enhancement in a-elements. 

(ii) Only for galaxies with H/3>2, unless their enhance- 
ment factor F and abundance ratios [Xei/Fe] (like [Ti/Fe]) 
are larger than the above limits, the presence of secondary 
star forming activity ought to be invoked. 

(iii) Looking at the position of models of constant 
metallicity and age but different F and/or [Xei/Fe] (for in- 
stance [Ti/Fe]), they scatter along a nearly verticalline. This 
implies that the metallicity relationship invoked for the bulk 
old population is no longer required. The vertical distribu- 
tion of the data is simply caused by the compensatory effect 
of different combinations of Z, F and [X^i/Fe] ([Ti/Fe] in 
our case). 

(iv) Secondary episodes of star formation are no longer 
a common feature to all galaxies, but an exceptional event 
limited to a small numb er of them. This ag r ees w ith the age 
distribution obtained bv iTantalo' TchioSdiooi). 

The a-Model is confirmed by the MonteCarlo simula- 
tions shown in the left and right panels of Fig. |7| In the left 
panel we show the case of old galaxies with no secondary 
activity: the bulk population spans the age range given by 
Ti = 13Gyr and r2 = 10Gyr with /3i = l, whereas the metal- 
licity and a-enhancement span the whole range for the pa- 
rameters Z, F, and [Xe;/Fe] ([Ti/Fe] taken as a measure of 
the natural width). Since there is no further star formation 
activity T3=10 Gyr and /32=0. The models essentially match 
the bulk of data, i.e. galaxies with H^^2, and yield a distri- 
bution in the H/3 vs. [MgFe] plane which is nearly vertical (no 
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Table 4. Relative errors for the theoretical simulations. 



(AH^/H^) 


(AMgb/Mgb) 


(AMga/Mga) 


(AC24668/C24668> 


(ANaD/NaD) 


(A(Fc>/(Fc» 


(A[MgFc]/[MgFc]) 


0.136 


0.059 


0.026 


0.009 


0.054 


0.073 


0.047 




Figure 7. Observational data and MonteCarlo simulations for the a-Model . Left Panel: The normal galaxies indicated 
by the open triangles and the shell- and pair-galaxies of iLonghetti et alJ iioOO) indicated by the filled squares. Central Panel; old 
galaxies formed in the time interval Ti=13 to T2=10 Gyr and no later star formation activity (/3i=l and /92=0). Each model galaxy has a 
metallicity and enhancement factor F randomly chosen from the whole range, i.e. 0.0085CZ^0.07 and O^FsCO.35. No universal enrichment 
law is considered. Thanks only to the different combinations of Z and F (different chemical compositions in general), the model galaxies 
span a large range of the observational values for and [MgFe]. Right panel: the same as before, but now a small percentage of a 
young population has been added, i.e. /3i=0.98 and /32 =0.02, whose age randomly varies in the interval O.l^TsSClO Gyr. The secondary 
burst of star formation has been included to get full coverage of the observational range in particular for those galaxies with H^>2.5. 
In both panels the meaning of the various lines (solid, dashed and dotted) is the same as in Fig. 1^ 



memory of the SSPs path). In the right panel we show the 
same but allowing for recent burst to occur, i.e. Ta^O.l Gyr. 
The burst intensity is for /32=0.02. But for the few galaxies 
clearly caught in the burst mode (those with H^>2), the 
two theoretical distributions are nearly identical and both 
fairly well reproduce the observational data. 

These simulations open the gate to an interesting al- 
ternative explanation, i.e. that the large scatter in and 
[MgFe] is predominately caused by a spread in the chemical 
parameters metallicity Z and enhancement factor F rather 
than metallicity and age in the bulk population of a galaxy. 
Secondary activity of star formation is unavoidable only for 
a minority of objects. 

The scatter in Z, T and also individual [X^i /Fe] of the 
dominant old stellar component could be attributed to dif- 
ferent kinds of star formation at the very early epochs, per- 
haps related to the physical conditions in the proto-galaxy 
affecting not only the intensity and duration of the star for- 
mation process, but also the initial mass function of the stars 
and the abundance ratios in turn. We will touch upon this 
point later. 



6 DISCUSSION 

Testing the a-Model. To corroborate the suggestion that 
the scatter in the vs. [MgFe] plane for most galaxies is 
due the chemical parameters F, Z, and perhaps [Ti/Fe] we 
perform here a few ad hoc experiments. 

First of all, we test a model with a very narrow age range 
for the old component and no subsequent star formation. 
The test is meant to isolate the sole effect of the chemical 
abundances. The model is characterized by the parameters: 
Ti^lSGyr, r2=12Gyr, /3i = l, r3=12Gyr, and /32=0. While 
no ULME is considered, the metallicities, Fs and [Xej/Fe]s 
(at present [Ti/Fe]) are let span the whole range. The sim- 
ulations are shown in left panel of Fig. |S] correlating Mgb 
with [MgFe]. The choice of these two indices is based on 
the notion that [MgFe] is nearly independent on F, whereas 
Mgb does. Both are sensitive to the age and metallicity even 
if some degree of degeneracy is present. On purposes we 
avoided H/3 because of its equal sensitivity to all the param- 
eters in question. In this plane we also displays the theoret- 
ical areas corresponding to the different combinations of F 
and [Ti/Fe], metallicities, and ages. Each area is bounded by 
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Figure 8. MonteCarlo simulations of the type shown in Fig. U] 
but with a narrower age range for the old component. They are 
aimed to better constrain the metallicity and enhancement factor 
suited to the data. Left Panel: models with only the very old 
component. Right Panel; models with both old and young stel- 
lar populations. In both panels the simulations are indicated by 
open circles and all other symbols have the same meaning as in 
Figs.|51and[3 In this pi ane real galaxies seem to fall in the ranges 
0^r<0.3 and 0.008<Z<0.05. 
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Figure 9. The same as in Fig. |S]but for the (Fe) vs. Mg2 plane. 
Now the data seem to be compatible with 0.35^r^0.50 and 
metallicities extending up to Z=0.07. 



the SSPs with the lowest (Z=0.008) and highest (Z=0.070) 
metallicities, the two heavy lines along which four values of 
the age are marked (2, 3, 8, and 13Gyr). It is worth notic- 
ing the good resolving power of Mgb only for F and large 
insensitivity to all remaining parameters. Different values 
of [Xej/Fe] ([Ti/Fe] in this case) have in practice no effect. 
Furthermore, the age and metallicity are nearly degenerate. 
The open c ircles are the sim ulations, the filled squares are 
the data bv lLonghetlj_etj1 j|200Q,') . and finally the open tri- 
angles are those bv lTraeeil ^ll997^ ■ Secondly, we re-examine 
in the same plane the simulations shown in Fig. |7| in which 
a wider age rage for the old component is adopted. All re- 
maining ingredients are the same as above. These model are 
shown in the right panel of Fig.|S| Comparing data with the- 
ory, we would conclude that in both cases F in the range to 
0.35 and Z ^0.07, likely 0.05, are best suited, the only dif- 
ference between the two type of models is that in presence of 
a burst the left-lower corner of the panels is populated (i.e. 
Mgb ^2 and [MgFe] ^3). In any case the large dispersion of 
the bulk population in the vs. [MgFe] plane is once more 
compatible with all galaxies being old and spanning a large 
range of Fs and Zs. The large age dispersion (say from 13 to 
5 Gyr) and the occurrence of a later burst of activity are no 
longer unavoidable results and/or hypotheses. Only a few 
galaxies seem to require a later burst of star formation. 

However, the upper limits suggested by the Mgb vs. 
[MgFe] plane are not firmly established as different con- 
clusions would be derived by looking at another diagnostic 
plane such as (Fe) vs. Mg2 shown in Fig. |5| All the symbols 
have the same meaning as in Fig. |H| From this diagram one 
would indeed conclude that F falls in the range 0.35 to 0.5, 
that high metallicities (up to Z—0.07) are likely to occur, 
and finally that galaxies with F^O and metallicities larger 
than solar or so are not there. Since Mg2 seems not to de- 
pend on F and (Fe) is likely more sensitive to both F and 
Z than other indices, we would favor the hints arising from 
this diagnostic plane with respect to the previous ones. 
Two Template Galaxies. It may be worth of interest here 
to closely inspect some of the galaxies in the samples paying 
major attention to those with unusually high H/j and/or 
very low [MgF e ]. To this aim th r ee sou rces have been used: 
iTraeeJ iiggTl) ; iLonghetti et alJ ll200(fl for the i ndices and 
central velocitv dispersions. iBursteh^taD I^Q^) for colors, 
masses, luminosities. The data are listed in Table|S] We have 
not considered RR24a and RR24b because of their negative 

Without going into a detailed discussion of all these 
objects, we may still make several points: 

(i) Adopting (B-V) = 1.00 ±0.05 as the typica l mean 
value for early- type galaxies iBurstein et al.lll997^ . those 
with strong tend to be bluer than normal galaxies by 
A(B-V)=0.05 to 0.10 up to the extreme case of NGC 3156 
with 0.28. Furthermore there seems to be a correlation be- 
tween (B-V) and H^: (B-V) gets redder at decreasing H/j. 

(ii) Looking at the time dependence of (B-V) (other 
blue colors have the same behavior) shown in Fig. blowing 
to the very fast decline past the occurrence of the burst, 
colors in the range 0.8 to 0.9 are recovered after 1.5 Gyr or 
so (the exact value depends on the intensity of the burst), 
which is the age of the young component 

(iii) From the position of a galaxy suspected to be in 
the burst or post burst activity because of its unusual in- 
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Figure 10. Age estimate of the last episode of star formation 
in a few galaxies with unusual values of and/or [MgFe] (Left 
panel) and Mg2 (Right Panel). The galaxies indicated with differ- 
ent symbols are listed in the Right Panel. In each panel we plot 
three SSPs with different metallicity, Z=0.008, 0.019, and 0.040 
as indicated, along which five values of the age are marked, i.e. 
13, 8, 5, 3 and 1 Gyr, and the corresponding lines of constant age 
are drawn (thin dotted lines) . Finally, we show the locus expected 
for composite galaxies, in which to a population with the age of 
13 Gyr and metallicity Z=0.019, an burst of arbitrary age and dif- 
ferent intensity is added. The metallicity of the young component 
is however the same as the old one (Z=0.019). The intensities are 
/32=0.005 (long dashed line), 0.01 (heavy dotted hne), and 0.02 
(solid line). Along each curve the age increases from the minimum 
plotted value 0.1 Gyr to 13 Gyr. The position of M32 is compatible 
with various combinations of Tb and /3, for instance Tb —1 Gyr 
and ,32=0.005 or Tg ~2-3 Gyr and /32 in the range 0.01 to 0.02. 
At given /32, the burst age will increase if higher metallicities for 
the young component are considered. 

dices in planes for composite systems like those presented in 
Fig-ISlon one hand we may get another estimate of the burst 
age, on the other hand we may check the consistency of the 
whole picture. For the sake of illustration, we discuss here 
the case of NGC 5018, a giant elliptical galaxy with total 
mass M=1.24xlO"M0, luminosity Lb = 5.43 x 10^°Lq, 
central velocity disper sion So=223 km/s. H/3=2.30, and 
Mg2=0.211. Long ago. iBertola et ahl il993l) have classified 
NGC 5018 as a giant, chemically unevolved galaxy consisting 
of M32-like stars, and^eonardi & Worthcy (2000.) have sug- 
gested that the central regions of NGC 5018 are dominated 
by an intermediate-age p opulation. Along th e same vein, is 
the very recent study bv lBuson et al.l (|2004| . private com- 
munication) who have stressed the close similarity between 
M32 (NGC 221) and NGC 5018. The similarity concerns the 
indices as shown by the entries of Table |51 the (B-V) colors , 
and the spectrum from 2000A to 4700A lBuson et al.l20ol) . 
They have also made the point that since M32 (at least in 
the central regions) is commonly considered to be dominated 



by a population whose age is about 3 Gyr (lFreedmanlll989l. 
1X9921: iTraeer et al.ll200o'bly: iLonghetti et al.ll2000h the same 
should occur for NGC 5018. However, the reality may be 
more articulated than this simple scheme. Studies of the 
stellar content and other properties of M32 have pointed 
out that the bulk population could be as old as that of 
globul ar clusters, e.g. see the discu ssion bv [ Bre ssaii et all 
ll99^ and ref erences therein, and bv lGrillniai^^Lri^9^r ^ 
iBressan et all |l994 , 1996) have argued that the young stel- 
lar component in M32 could result from a recent burst of 
activity superposed to the pre-existing, old stellar content. 
Is NGC 5018 experiencing the same? To check this possibil- 
ity, in Fig.Iiniwe plot M32, NGC 5018 and other galaxies of 
Table 13 on to the H/j vs. [MgFe] and Mg2 planes and com- 
pare them with single SSPs of difi'erent age and metallicities 
and composite SSPs in which bursts of different intensity 
and age are superposed to an old dominant component. In 
the first case we recover the standard solution, i.e. M32 is 
an object whose stellar content is about 3-4 Gyr, whereas in 
the second case M32 is an old object (say 10-13 Gyr) on the 
top of which a young component (1-2% of the total mass) is 
present. The age of this is about 1 Gyr. An old galaxy, which 
underwent a burst of activity say 1 Gyr ago, would appear 
as a brand newly formed galaxy with the age of about 3- 
4 Gyr. The same analysis can be extended to other galaxies 
of Table |3 A special remark is due to the three galaxies 
with very low [MgFe]: we exclude RR317b because it shows 
evi dence of emission [OII ](3712A) as already pointed out 
bv ILonghetti et al.l fcOOOh . The case of NGC 3156 is triv- 
ial because of its strong H/3: it is indeed fully compatible 
with theory if undergoing a burst of star formation. We 
are left with RR62a for which not very much can be said: 
it follows the same ^e-Re relationship of ordinary galax- 
ies, but it deviates from the Hamabe-Kormendy relationship 
jHamabe fc Kormendvl ll987'l. It could be a old, low metal- 
licity galaxy which according to lGorgas et alJ il993l) should 
fall into the bottom left corner of the H^g vs. [MgFe] plane. 
Alternatively it could be old galaxy which is just starting 
(or undergoing) a burst of star formation as indicated by 
the results for composite SSPs shown in Fig. |3 
The Coma EGs. We have already mentioned that EGs in 
Coma Jj0rgensenlll999D have systematically higher values of 
Hp than the local galaxies. The minimum values is Ha=1.5. 
Repeating the same analysis we have made for the iTrageJ 
|l992J and Longhetti et aL (^OOO,) samples we would come 
up with the conclusions: (i) The a-Model may hold good 
even in this case, (ii) The Coma galaxies have a higher degree 
of enhancement than the local ones. A rough estimate is 
r ^0.3. Suppor ting this idea is the es timate for the age of 
Coma's EGs bv lPoggianti et all i200ll) who find a mean age 
older than about 9 Gyr; the e s timat e of the ratio [Mg/Fe] by 
|j0rgensenl lll997f ): Ijorgenseni lll999l) who finds [Mg/Fe]=0.3 
to 0.4 as the central velocity dispersi on log Sq increases o f 
0.4 dex; and finally the recent study bv lMehlert et alJ (|2003h 
who find [a/Fe] ratios in the range 0.15 to 0.4. The mean 
degree of enhancement in EGs of the Coma Cluster seems 
to be at least 0.1 dex higher than in the local field galaxies. 
Q-Enhancement in EGs: Data and Theory. The oc- 
currence of the a-enhancement in early-type galaxies is cur- 
rently attributed to the duration of the star forming period 
and the different contribution to che mical enrichmen t by 
Type la and Type II supernova (see iMatteuccil 119971 . for 
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Table 5. Selected galaxies: seven indices on the Lick system; M is the mass in solar units inside the effective radius; Lg is the B-luminosity 
in solar units; So is the central velocity dispersion in km/s. 
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a review of the subject). In brief, a-elements are mainly 
ejected by Type II together with some Fe, whereas Fe is 
essentially expelled by Type la supernova (via the car- 
bon ignition in binary White Dwarfs reaching the Chan- 
drasekhar limit by mass accretion). As long as Type la 
supernovae do not intervene in a significant manner, the 
chemical composition of the gas and newly formed stars 
in turn will be enhanced in a-elements. Since Type la 
supernovae are expected to star t contaminating the inter - 
galactic medium than Type II jGreggio fc Renzinil Hoil . 
in the framework of the standard supernova driven galac- 
tic wind model by iLarsonl (^T^) and the standard initial 
mass function, the time scale of star formation and galactic 
wind must be shorter than about 0.5 Gyr not to decrease 
the initial [a/Fe]>0 (when a-elements are mostly pro- 
duced) to [a/Fe\^0 (when iron is predominantly ejected). 
In other words, to reproduce the observed trend of the 
[a/Fe]-mass relationship, the total duration of the star 
forming activity ought to scale with the galaxy mass ac- 
cording to AtsF OC . This is the main drawback of the 
standard model because it requires a mass-star-formation- 
timescale correlation which is the opp osite of what implied 
by the Color-Magnitude relationship (Bower ot al. 'l992a') 
as amply discussed inlB rcssan ct al. (1996); Chiosi (2000); 
IChiosi fc Garrard J200a and lTantalo fc Chiosil ll20Q2i) . This 
point of difficu lty has been overc o me by the N-Body Tree- 
SPH model of IChiosi fc Garrard tOO'il . In brief (i) Inde- 
pendently of the total mass, galaxies of high initial density 



undergo a prominent initial episode of star formation fol- 
lowed by quiescence, (ii) The same applies to high mass 
galaxies of low initial density whereas the low mass ones un- 
dergo a series of burst-like episodes that may stretch over 
a considerable fraction of their lifetime, (iii) The mean and 
maximum metallicity increase with the galaxy mass. There- 
fore these models can account for the GMR of EGs. (iv) 
Finally, th e occurrence o f galactic winds does not follow 
the simple iLarsonl il974ft model, but takes place continu- 
ously in lumps of gas as soon as their thermal kinetic en- 
ergy exceeds the gravitational e nergy of the ga laxy. In other 
words the scheme on which the iLarsoiil il974) model rests, 
i.e. - longer star formation period - higher metallicity - 
lower degree of a-enhancement at increasing galaxy mass, 
is reversed but for the metallicity. Therefore the GMR and 
Q-enhancement constrai nts ar e met at the same time. See 
iGhiosi fc Garrard tOO'j) and iTantalo fc Ghiosil J2002h for 
all details. Given these premises, the existence of various 
degrees of Q-enhancement in different EGs finds a natural 
explanations thanks to the sensitivity of the star forming 
process and its duration to the initial conditions (mean den- 
sity) and total mass of the proto-galaxy. By the same token 
we may also explain the systematic higher abundance ratios 
[a/Fe] passing from field to cluster galaxies. 
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7 WHAT IS THE MAXIMUM COLOR 
DISPERSION OF EC'S? 

Our aim here is to check at what extent a past episode of 
star formation would reflect on to broad-band colors, such 
as (B-V), of the host galaxy as we see it today"^. To this 
aim we adapt the toy model of galaxy formation we have 
been using in so far. The analysis below applies both to the 
monolithic (isolation) and hierarchical scheme. 

For the sake of simplicity let adopt a model with the 
minimum number of parameters: solar scaled abundances, 
fixed age for the old stellar component, here indicated as 
Tf —12 Gyr, flxed mean metallicity. With no further stellar 
activity, passive evolution of the galaxy would imply that at 
the present time (B-V)~1.00±0.05 whereby the uncertainty 
reflects the spread in metallicity and abundance ratios. 

Now suppose that at a certain age Tb the model galaxy 
undergoes an additional episode of star formation of short 
duration and arbitrary intensity P2 (this immediately fixes 
also the intensity /3i of the old component, their sum being 
equal to unity). The burst can be of internal origin (e.g. gas 
left over by the previous activity and then turned into stars) 
or the result of a merger with some gas-rich object. In the 
first case, the galaxy will be made by an old and a young 
component. In the second case, the galaxy will consist of two 
sub-units harboring old populations (in principle of different 
age and metal content) plus the younger component born 
during the burst. 

The color evolution of the composite galaxy is followed 
up to the present time and its broad-band colors are tested 
against observations. The results are displayed in the large 
panel of Fig. 1111 as a function of the age and intensity of 
the secondary activity. Let us now read off the age Tr, at 
which a galaxy which underwent a burst of star formation of 
intensity 132 at the age Tb, is able to recover the same color 
it had before the burst. This is simply given by the inter- 
section of any line labelled (32 with the horizontal bar. The 
ages Tr as function of the /?2 are shown in the small panel of 
Fig. 111! It turns out that for many combinations of Tb and 
(32 , the resulting (B-V) color would be too blue as compared 
to the typical colors of EC's, (B-V) = 1.00±0.05. Stellar ac- 
tivities engaging 5-10% of the total mass and taking place 
as early as 5-6 Gyr ago would be detectable. The situation 
gets even worse for higher /32 and/or lower Tb- This implies 
that only remote or minute star forming events are allowed. 
Therefore, either mergers occur very early on in a galaxy's 
history or only captures of small bodies and little compan- 
ion star formation can take place at later epochs. This is a 
rather strong constraint that should be taken into account 
by any model of hierarchical assembling of big galaxies. It 
is plausible to suggest that the main body of a galaxy was 
assembled early on and that all subsequent activity we in- 
fer from the observational data is limited to the capture of 
small objects (likely satellites of the dominant galaxy). 



^ P reliminary studies of the same subject are by IChiosil I2OOCII) 
and lChiosi fc Carrarol i2002h 
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Figure 11. Simulations of the (B-V) color of EGs in which a 
burst of star formation of arbitrary intensity ad age Tg is su- 
perposed to an old population with the age of 13 Gyr. In these 
simulations for both components we assume Z=0.019 and r=0. 
The burst intensity is /32=0.001, 0.002, 0.005, 0.01, 0.02, 0.08, 
0.10, 0.20, 0.30, 0.40, and 0.50, going from bottom-left to top- 
right as indicated. The horizontal bar visualizes the typical color 
of EGs with the uncertainty A(S - V) = ±0.055. The insert 
shows the time Tr in Gyr elapsed from the epoch of the burst to 
the stage at which the composite galaxy recovers the typical red 
colors of EGs as a function of the burst intensity. Tr is given by 
the intersection of each color curve with the horizontal bar. 

8 SUMMARY AND CONCLUDING REMARKS 

The present-day challenge with EGs is to unravel their for- 
mation and evolution history. In a simplified picture of the 
issue, the problem can be cast as follows. Do EGs form by 
hierarchical merging of pre-existing sub-structures (maybe 
disc galaxies) made of stars and gas? Was each merging 
event accompanied by strong star formation? Or conversely, 
they originate from the early aggregation of lumps of gas 
turned into stars in the remote past via a burst-like episode 
ever since followed by quiescence so as to mimic a sort of 
monolithic process? Even if the two alternatives seem do 
oppose each other, actually they may concur to shaping the 
final properties of EGs as seen today. The question is to un- 
ravel the signature of the forming mechanism from the obser- 
vational data. To this aim we have examined the line absorg- 
tion in dices on the Lick system of normal, field EGs of Trageg 
l|l997|) and of interact ing EGs (pair- and shell-objects) of 
iLonghettT et al.l ll2000D . Occasioi ially we have also looked at 
the EGs in the Coma Clusters of lJorgensenl (^9^ and other 
more specific cases taken from literature. The bottom line 
of this study was (i) to check whether normal (quiescent) 
and interacting have a different behavior in the popular di- 
agnostic planes such as vs. [MgFe] (and others); (ii) to 
seek whether the signature of interaction may mirror in some 
specific changes of the indices that could un-equivocally hint 
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for bursts of stellar activity; (iii) to evaluate the intensity of 
those bursts or secondary episodes of star formation; (iv) 
to explore whether other alternatives can exist, i.e. distinct 
from obvious ones resting on large age range and/or bursts 
of star formation (from one to several) at various epochs. 

From the various observational issues we have been ex- 
amining in so far, we gather the following picture: 

(1) Both normal, field and interacting galaxies have the 
same scattered distribution in the Hp vs. [MgFc] diagnostic 
plane even the interacting ones show a more pronounced tail 
toward high Hp values. This may suggest that a common 
physical cause is at origin of their distribution. 

(2) The distribution of normal and interacting galax- 
ies is smoothly elongated in H/3. While for the interacting 
objects an easy, however not fully satisfactory, explanation 
can be found invoking a late burst of activity, in the case of 
normal galaxies an explanation can be found only invoking 
a large age range for the bulk population together with the 
existence of an universal law of metal enrichment. Adding 
a late burst of star formation to an old system without in- 
voking the large age range and the universal law of chemical 
enrichment would not yield the desired trend. Therefore not 
only normal galaxies have to be build up over a large time 
interval but also at increasing total metallicity. The occur- 
rence of late star forming episode is a sort of optional. 

(3) More specifically, a typical model based on clas- 
sical SSPs with solar partition of elements and where the 
burst is superimposed to an old and coeval population is 
not able to reproduce the smooth distribution of galax;ies in 
the H/3 vs. [MgFe] plane. This kind of model would predict 
an outstanding clump at low values, contrary to what is 
observed. Models in which star formation lasted for a sig- 
nificant fraction of the Hubble time (4Gyr ^toid^lG Gyi) 
better match the observed diagram. In this context, the pe- 
culiar, almost vertical distribution of galaxies (normal, shell- 
and pair-objects) in the Ha vs. [MgFe] plane is interpreted 
as the trace of the increase of the average metallicity accom- 
panying all star forming events. This could be the signature 
of a metal enrichment happening on a cosmic scale. 

(4) However, the above scheme is firstly too demand- 
ing because of the many ad hoc ingredients that have to 
introduced, secondly it neglects important effects given by 
the observationally grounded hint that the stellar content of 
EGs is likely enhanced in a-elements with F ranging from 
0.1 to 0.4 dex. We would like to propose a new scheme, in 
which the bulk dispersion of galaxies in the H^ vs. [MgFe] 
plane is caused by a different mean degree of enhancement. 
Indeed two old galax;ies of the same age (say 13 Gyr) but dif- 
ferent r and [Xej/Fc] for some important species (C,N,0... 
Ti) would have different values of H/3. For instance passing 
from a 13 Gyr old galaxy with r=0 to an object with the 
same age but r=0.3 would increase H^ by as much as 0.2- 
0.3 or more (depending on the abundance ratios for some 
specific elements like Ti). The effect is comparable to the 
mean observational dispersion. The majority of EGs can be 
pretty old (10-13 Gyr). Furthermore, neither large age range 
nor universal enrichment law axe required. The nearly ver- 
tical smooth distribution for coeval galaxies is secured by 
the compensatory effect of different combinations of Z and 
P. Finally, the occurrence of a late burst of activity is not 
an un-avoidable ingredient of the recipe, but a rare event 
interesting only those galaxies with very high Hp (roughly 



> 2.5). The possibility that EGs span large ranges of F and 
mctallicities but narrow ranges of ages for the bulk popula- 
tion favors the monolithic scheme and can be derived from 
N-Body Tree-SPH models of galaxy formation as a conse- 
quence of the dependence of the star formation efficiency 
and temporal history on the initial mean density and total 
mass of the proto-galaxy. 

(5) As far as we can tell, from simulations of broad- 
band colors, EGs are compatible with the notion that the 
bulk stars have formed in the remote past. Galaxy mergers 
and companion star formation in a recent past are not likely, 
unless the intensity of the secondary activity is very small 
(i.e. engaging less than a few percent of the total mass). 
Were merging the only possible mechanism to form massive 
EGs, this should have occurred in a remote past (half of 
the Hubble time at least). In any case merging of smaller 
units without star formation even in a recent past cannot 
be excluded. Is it likely? 

(6) Finally, prolonged or secondary stellar activity seem 
to be also more probable in field and loose groups EGs 
than in those belonging to compact groups and clusters. The 
physical cause could once more be the mean density of the 
environment out of which proto-galaxies are formed. 
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